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Magnetic field tuning and quantum interference in a Cooper pair splitter 
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Cooper pair splitting (CPS) is a process in which the electrons of naturally occurring spin-singlet 
pairs in a superconductor are spatially separated using two quantum dots. Here we investigate 
the evolution of the conductance correlations in an InAs CPS device in the presence of an external 
magnetic field. In our experiments the gate dependence of the signal that depends on both quantum 
dots continuously evolves from a slightly asymmetric Lorentzian to a strongly asymmetric Fano- 
type resonance with increasing field. These experiments can be understood in a simple three - site 
model, which shows that the nonlocal CPS leads to symmetric line shapes, while the local transport 
processes can exhibit an asymmetric shape due to quantum interference. These findings demonstrate 
that the electrons from a Cooper pair splitter can propagate coherently after their emission from 
the superconductor and how a magnetic field can be used to optimize the performance of a CPS 
device. In addition, the model calculations suggest that the estimate of the CPS efficiency in the 
experiments is a lower bound for the actual efficiency. 

PACS numbers: 74.45.+c, 73.23.-b, 73.63.Nm, 03.67.Bg 


In the Cooper pair splitting (CPS) process the elec¬ 
trons of the Cooper pairs in a superconductor are sep¬ 
arated spatially using two quantum dots (QDs) coupled 
in parallel to a central superconductor contact (S) in a 
three-terminal geometry m, see Fig. la. The Coulomb 
repulsion on the QDs and the quasiparticle energy gap of 
the superconductor enforce the electrons to separate into 
different normal metal electrodes (N1 and N2). Since 
Cooper pairs are spin singlet states, such devices could 
serve as a source of nonlocal spin entangled electron pairs. 
Similar geometries are also relevant in the search for Ma- 
jorana bound states [4], both, in a local S-N junction ex¬ 
periment or in a three-terminal setup, where an increase 
in CPS efficiency might serve as a signature of the elusive 
exotic states [5]. 

In a series of recent experiments on semiconducting 
nanowires (NWs) [6] 9 . carbon nanotubes [TQ14T2] . and 
graphene m, CPS was demonstrated by positive con¬ 
ductance correlations between the currents from S into 
the two N-contacts. In these experiments, external mag¬ 
netic fields were solely used to suppress the superconduc¬ 
tivity for control experiments, but not as a parameter 
to tune CPS. In addition, most experiments were inter¬ 
preted in terms of an incoherent picture with independent 
transport mechanisms, only coupled by the quantum dot 
(QD) dynamics 9. ITT]. 

Here we report experiments in a NW-based Cooper 
pair splitter with a Nb superconducting electrode. The 
large critical magnetic fields of Nb allow us to explore 


CPS up to ~ 1 T. We find that the conductance correla¬ 
tions cannot only manifest as symmetric peaks and dips 
in the gate dependence, but also as strongly asymmet¬ 
ric shapes reminiscent of Fano resonances. We interpret 
the experimental results in a minimal model that incor¬ 
porates the superconducting proximity effect, the tun¬ 
nel coupling between the QDs and quantum interference 
(see Fig. lb). Interference results in asymmetric Fano- 
type features for local transport processes, whereas CPS 
produces Lorentzian shaped contributions to the conduc¬ 
tances. These effects can be strong enough to obscure 
the positive correlations due to CPS. The interference 
features in our experiments suggest that electrons can 
be transmitted coherently through a Cooper pair splitter 
device, which is a fundamental prerequisite for testing 
Bell’s inequality in a beam mixer setup HU. Since of- 



FIG. 1: (Color online) (a) and (b) Schematics of a Cooper pair 
splitter and the corresponding 3-site model. The dashed and 
solid gray arrows represent two different interfering single¬ 
electron paths from S to N2, which can result in Fano-type 
conductance resonances. 
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ten the spin-coherence time is longer than the spatial co¬ 
herence, our results are also encouraging for the various 
propositions to detect spin entanglement H5HEJ. 

Figure shows a schematic of a CPS device. Us¬ 
ing electron-beam lithography, an array of twelve 40 nm 
wide local bottom gates (4/18 nm Ti/Pt, ^ 60 nm spac¬ 
ing) were fabricated on a silicon substrate, electrically 
insulated by ^ 25 nm SiN x . A single NW (~ 70 nm di¬ 
ameter) was then deposited perpendicular to the gates 
using micromanipulators. The NWs were grown by solid- 
source molecular beam epitaxy [20], using an optimized 
process to suppress stacking faults [21]. A 330 nm wide 
and 110 nm thick superconducting Nb contact in the cen¬ 
ter and two normal metal electrodes (7/95 nm Ti/Au) at 
the ends of the NW were fabricated consecutively by con¬ 
ventional electron beam lithography, using an ammonium 
sulfide passivation [22] to remove the native NW oxide. 

The experiments were carried out in a dilution refrig¬ 
erator at a base temperature of ^ 50 mK. The QDs on ei¬ 
ther side of S were each induced in the NW by a negative 
voltage applied to 2 neighboring bottom gates, of which 
one was also used to tune the chemical potential of the 
respective QD. We label the voltages on these two tuning 
gates as V g i and V g 2 and use the terms ’local gate’ and 
’far gate’ to distinguish the two gates when discussing a 
specific QD. A sinusoidal voltage U( ac ) « 10 /xV applied 
to S results in the simultaneously recorded currents l[ ac ^ 
and I 2 in the contacts N1 and N2, which were held at 
carefully nulled potentials. We define the differential con¬ 
ductance through each QD as G{ = /U c )/y( ac ), both of 
which show well-defined, uncorrelated Coulomb blockade 
diamonds [23] . 

Figures 2a shows G\ and G 2 as a function of the two 
QD tuning gates for B = 0. The resonance amplitudes 
of both QDs are increased when both are on resonance 
(resonance crossing), which results in a positive corre¬ 
lation between the conductance variations. In contrast, 
in the measurements at B = IT in Fig. 2b, the con¬ 
ductance maxima do not occur exactly at the expected 
resonance crossings, but vary in position, which results in 
a conductance not symmetric with respect to a crossing. 
Figs. 2c and d show cross sections taken on top of the 
a QD resonance as a function of the respective far gate 
for a series of magnetic fields. Away from a resonance 
crossing the resonance positions and amplitudes are in¬ 
dependent of the respective far gate, while near a crossing 
both vary with gate voltage and exhibit different symme¬ 
tries with respect to the resonance crossing (dashed lines) 
with increasing field. For B = 0 both resonance mod¬ 
ulations are dominated by a roughly Lorentzian shape 
with a small minimum at the low-voltage side of the res¬ 
onance crossing. As we increase the magnetic field, G\ 
and G 2 evolve differently: while G\ changes only little up 
to B « 400 mT and then becomes broadened at higher 
fields, the variation of G 2 first evolves into a dip/peak 
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FIG. 2: (Color online) (a) and (b) G 1 and G 2 as a function 
of Vgi and V g 2 for B — 0 and B — IT, respectively, (c) and 
(d) Gi and G 2 as a function of the respective far gate, V g 2 
and V g i, for a series of magnetic fields B with the local gates 
set to a QD resonance, as indicated by black arrows in Fig. a 
and b. The curves are shifted vertically for clarity. 


structure, into a dip at intermediate fields and a max¬ 
imum at the low-voltage side of the resonance crossing 
for the highest fields. The finite-field curves are reminis¬ 
cent of Fano-type resonances, rather than Lorentzians. 
Directly comparing the B = 0 to the B = IT curve, 
one finds that the conductance maximum is shifted from 
the right to the left side of a crossing. We note that for 
most curves the maxima and minima in the conductance 
variations do not coincide with the resonance crossing. 

The resemblance of the experimental conductance vari¬ 
ations to Fano-type resonances suggests that at least one 
of the transport mechanism is prone to quantum interfer¬ 
ence, which cannot be described in the incoherent mod¬ 
els of previous works 0 m- To account for the ob¬ 
served interference patterns also requires to go beyond 
earlier coherent two-dot models uni m Eg, where the 
QDs are coupled by direct tunneling. Here we introduce 
a 3-site model, shown schematically in Fig. [TJd, in which 
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each QD is represented by a single spin-degenerate level 
with Coulomb interactions and coupled coherently to a 
central NW segment below S and to the leads N1 and 
N2. The central segment is also modeled as a single level 
coherently coupled to S. The Hamiltonian (without the 
coupling to the leads ESI) reads 

Hq — €i f crTli f< j T T ^ ^ T h.C., 

a,i cr,i^m 

where i = l,ra,2 label the left, middle and right levels 
with the on-site energies = e^(F>) + cr^H/2, com¬ 
prising orbital and Zeeman shifts. For the g-factors 
gi we use typical values in the range of \g\ = 5 — 15 
[261127]. ti m corresponds to the hopping amplitudes from 
the QDs to the central site. To simulate the supercon¬ 
ducting proximity effect we set the Coulomb interactions 
on the central region to zero, U m = 0, which is justified 
by the large size and the screening by the superconduc¬ 
tor. This assumption allows the occupation of the cen¬ 
tral level by two electrons, e.g. a Cooper pair from S. 
We assume that only this central region is coupled to 
S. 7 £, with i = l,m,2, are the tunneling rates from the 
three sites to the respective lead, see Fig. lb. We cal¬ 
culate the linear transport characteristics using an equa¬ 
tion of motion approach for the electronic Green’s func¬ 
tions [23, 28] . The conductances can be decomposed as 
Gi = Gi OC ;i+Gcps where the local contributions comprise 
local Andreev reflection (LAR) and single quasiparticle 
tunneling, Gi oc; i = G L AR ; i + G qp;i . G C ps denotes the 
contribution due to CPS. 

The 3-site model qualitatively reproduces the experi¬ 
ments in Fig. 2. The corresponding calculated transport 
characteristics are plotted in Fig.[3]for typical experimen¬ 
tal parameters. Here, — q corresponds to the changes in 
gate voltage (up to a lever arm factor). The model ac¬ 
counts for the shifts of the resonance maxima and min¬ 
ima with increasing magnetic field, the relatively weak 
change in G\ at low fields, as well as the peak/dip/peak 
transitions in G 2 with increasing field. 

To intuitively understand the structure of the con¬ 
ductance variations and the field evolution, we plot in 
Fig. [S5} i the local and CPS contributions to the total 
conductance G 2 for B = 0 and B = 1.2 T. We find that 
Gioc of both QDs often exhibits a strong asymmetric dip 
as a function of the respective far gate voltage, while 
Gcps generally results in a Lorentzian peak, not neces¬ 
sarily at the same position. The physical reason for these 
characteristics is the impact of interference: for the local 
processes into N2, the lowest order non-zero contribution 
comprises two different paths. Either an electron tunnels 
directly from the central region to QD2 and to lead N2, 
or it reaches QD2 and N2 after an excursion to QD1 and 
tunneling through the middle region, see Fig. [TJd. Since 
the transmission phase acquired on a QD state can vary 
by a large fraction of 7 r [29], the two paths can interfere 
constructively or destructively, depending on the relative 



(b) 


Gi(2e 2 /h) 





-e 2 (A) -ei (A) 


1 

0.52 
0.5 Q 

to 

0.48 IT 

0 . 4 " 

0.44 


FIG. 3: (Color online) Conductances through QD1 and QD2 
in the three-site model with parameters chosen to qualita¬ 
tively reproduce the experiments in Fig. [ 2 ] (a) and (b) Gi 
and G 2 as a function of both QD level positions, (c) and (d) 
On-resonance conductance variations of QD1 and QD2, re¬ 
spectively, as a function of the level position of the other QD, 
for increasing values of magnetic field. The parameters used 
in the calculations are e m — —0.03A, 2£i m — — 0.1 A, 

U=3A, and the coupling to the leads 71 = 0.15A, 72 = 0.2A, 
and 7 m — 0.11 A. 


level positions, resulting in a minimum or a maximum 
in the local contribution and in an asymmetric conduc¬ 
tance variation. In contrast to the local currents, the 
lowest-order non-zero CPS contribution stems from the 
direct tunneling of one electron into each QD, which is 
not affected by interference [23] . The sum of the local 
and nonlocal processes can then result in the observed 
Fano-type line shapes. 

The inversion of the line shape symmetry with mag¬ 
netic field, shown in Fig. 4a, can be understood as follows: 
the orientation of a Fano resonance depends on a generic 
shape-parameter [3], which is given here essentially by 
e m [23]. When e m is below (above) the Fermi energy, the 
asymmetry of the local contribution has a maximum at 
a more positive (negative) gate voltage than the dip. In 
the experiments of Figs. and d we find a transition 
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FIG. 4: (Color online) (a) Total conductance through the 
QD2, G 2 (solid line), and the decomposition into the local 
(dashed line) and CPS (dotted line) contributions, as a func¬ 
tion of —ci for B — 0 (dark) and and B = 1.2T (green), (b) 
G 2 vs. V g i for different coupling strengths 7 m to the super¬ 
conductor. The curves are offset vertically for clarity. The 
value at which 7 m = ti m is pointed out by an arrow, (c) and 
(d) Correct and estimated CPS efficiencies, x and Xe, f° r the 
numerical results in (b). 


from a peak/dip to a dip/peak structure at B « 0.6 T, 
which we thus identify as the field at which the posi¬ 
tion of the central level is shifted from below to above 
the Fermi energy by orbital and Zeeman shifts. We note 
that an asymmetric shape also shifts the maxima of the 
total conductance with respect to the resonance cross¬ 
ing. We observe these characteristics in the experiments 
(Figs. 2c and d) as well as in our calculations (Figs. 3c 
and d), where the maxima and minima do not occur at 
the dashed lines that indicate the crossings. 

As can be inferred from Fig. 4a, the CPS contribu¬ 
tion does not change significantly with B (for low enough 
B), but the local contribution can be reduced strongly. 
The latter is due to the lifting of the spin-degeneracy 
by the Zeeman shift, which results in a reduction of 
LAR due to the opposite-spin electrons of a Cooper pair 
tunneling through a (partially) polarized level. This 
leads to an increase in the CPS efficiency defined as 
X = 2Gcps/(Gi + G2 ) 0111 ]. 

The efficiency is relevant for most prospective applica¬ 
tions of a CPS device HUHH. In the above expression for 
X one usually estimates Gcps by subtracting the conduc¬ 
tance sufficiently far away from the resonance crossing. 


We call this experimentally obtained efficiency Xe- It was 
shown for an incoherent model that this procedure un¬ 
derestimates the actual CPS efficiency [9]. However, this 
is not obvious for coherent transport. In Fig. 4a we show 
in an example that the local contribution is not constant, 
but exhibits a minimum at the resonance crossing, which 
also suggests that x > Xe- To further illustrate this effect, 
we have calculated G 2 vs. c\ (i.e. V g \) in the 3-site model 
for a series of tunnel coupling strengths between S and 
the middle level, see Fig. 4b. We find that the conduc¬ 
tance minimum for weak couplings evolves into a max¬ 
imum for 7 m > ti m (black arrow). The reason for this 
transition is an additional contribution by states forming 
in the superconducting gap as increases. In Fig. 4c, x 
is plotted for all curves in Fig. 4b and is almost identical 
for all 7 m . For comparison, Xe is plotted in Fig. 4d for 
the same conductances, which demonstrates that Xe sys¬ 
tematically underestimates x (we have omitted negative 
values) [31] . 

In summary, we report experiments on a Cooper pair 
splitting device that exhibits characteristics of quantum 
interference, which are tunable by electrical gates and ex¬ 
ternal magnetic fields. To interpret these data, we intro¬ 
duce a 3-site model that reproduces the experiments on a 
qualitative level and allows the decomposition of the con¬ 
ductances into local and CPS contributions. Specifically, 
we show that quantum interference is relevant in such 
structures, which is required for any application of CPS 
as a source of entangled electron pairs. In addition, we 
discuss the impact of various device parameters and find 
that the experimentally estimated CPS efficiency is quite 
generally a lower bound for the actual value also in the 
coherent transport regime. The presented model might 
also be used to gain a new vantage point for other fun¬ 
damental effects in superconducting proximity systems, 
such as fractional fermions or Majorana bound states. 
For example, a similar setup could be used to detect the 
inversion of the gap in a Majorana wire by the change of 
symmetry in the local Andreev conductance as a single 
level crosses the Fermi energy. 
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SUPPLEMENTARY MATERIAL 
QUANTUM DOT CHARACTERIZATION 

The Cooper pair splitter device was equipped with 12 bottom gates (see Fig. la in the main text). Starting the 
numbering on the left, we applied negative voltages to the 2nd and 3rd gate to form the barriers for a quantum dot 
(QD) in the left arm (QD1), and to the 9th and 10th gate to form QD2 in the right arm. Gate 2 was then used to 
tune QD1 and labeled gl in the main text. Similarly, gate 8 was used to tune QD2 and labeled g2. Representative 
conductance maps of QD1 and QD2 as a function of the bias applied to S and the respective tuning gate voltage are 
presented in Fig. SI. 
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FIG. S5: (Color online) (a) and (b) Representative conductance maps of QD1 and QD2, respectively, as the function of the bias 
and the respective tuning gate voltage. The charging energies are in the range of E c « 2 — 3 meV for QD1, and E c « 8 — 12 meV 
for QD2. The FWHM (full-widths-half-maximum) are ~ 130 — 200 geV for both QDs. 


MODEL HAMILTONIAN FOR INTERACTING ELECTRONS: SELF-CONSISTENT APPROACH 

In this section we provide a description of the self-consistent method used for the calculations in the main text |lj. 
Model and Hamiltonian — We start by defining the model Hamiltonian for the three tunnel-coupled sites in our 
3-site model (see Fig. lb of the main text), 

Hq — ^ ^ G,<7^,(7 T T ^ ^ timdi ^dm^cr T h.C.. (1) 

cr,i cr,i^m 

As described in the main text, QD1 and QD2 are coupled to the normal leads 1 and 2, respectively, each described 
as non-interacting Fermi liquids with a Hamiltonian of the form 

Hi = £k a ]k a ^kcr- (2) 

k a 


a !kcr are the annihilation (creation) operators in for electrons in lead l with momentum k and spin cr. ti m is the 
tunneling amplitude between the central site and site i. 

The central site is coupled to the superconducting lead S, for which we use the BCS-Hamiltonian 

Hs = 5^(^ kC kt C -kt + ^k c -kt c kt) (3) 

k,cr k 

where c£. creates a fermion with momentum k and spin a. 

The coupling to the leads is given by 

H T = Y h (4,<T a i,k,<T + h.C.) +y^t s (d^Cka + h.C.) 

k,cr k,cr 

2 = 1,2 


( 4 ) 
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Tunneling from a QD to the state k in the normal lead l is described by the amplitude £/. We assume that the 
k-dependence of the amplitudes ti and t s can be neglected (wide band limit). Thus, we can define the tunnel rates 
between site i and lead i as 7 $ = npit? (for i= 1 , m and 2 ), with pi the density of states per spin in lead i. 

Green’s function of the interacting system — In the presence of the Coulomb interaction, it is not possible to obtain 
the exact Green’s function (GF). We therefore decouple the corresponding equation of motion (EOM) at the level of 
two-body GFs. Under this approximation the one-body GF 

G{t, t') = iO(t' - WOT) 


coupled to the two-body GF by [T] 


w(t,t r ) = -mmMty})), 


which can be written in the basis 

V> f 0 ) = 


n t 
a it 


dll ^mt d\i d2l 


^(t) = d u n l,t d lnt nm ’l d mirim, f < 4 f n 24 d 2\.n 2 ^ 


( 5 ) 

(6) 

( 7 ) 

(8) 


The dot operators are used in the Heisenberg picture labeled as £ t’. The 3-sites model leads to a 12 x 12 GF, but in 
the absence of a magnetic field we can use the time reversal symmetry to express the Green’s function in a reduced 
6 x 6 matrix form. 

The time derivative of Eqs. © and © yields the equation of motion (EOM) with the Fourier transform 

{Cj-ho-t- f )G(w) = 1 + UW(lj), and (9) 


W{u) = {cj-h 0 - C/)“M 


■s + r 


GH, 


( 10 ) 


where we omit terms involving Kondo-like correlations. We have used the matrices ho = J 2 i j=i m r ( — l) 2+1 £i 
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( 11 ) 


( 12 ) 


(13) 


where the coupling to the normal contacts gives rise to the self-energies 7 *, and the coupling to the superconductor 
results in the self energies 


/ x _ u-irjs 

9 ” ( " ) - 7 ” y a 2 - (u, - i,.)2 

a / \ A 

^indyU) — 7 m - - - : - — • 

V A 2 - (u - ir)s) 2 


(14) 

(15) 




















The first term represents the self energy due to the quasiparticles in the superconductor, while the second term is the 
anomalous self energy that induces a superconducting gap on the central site. 

Substituting Eq. (10) into Eq. (|9|, we obtain 


GV) = (V - S - f ) 


-1 


(16) 


where 


ffo = A - M” 1 (l + U(u -ho- U)~ l A^ . 


(17) 


Equation (16) needs to be calculated self-consistently because the matrix A contains parameters (n ^ cr and (di^d^)) 
that depend on the resulting G{uS). 

Linear conductance — In the linear response regime, with symmetric biasing and T -A 0, the conductance through 
QD l contains four contributions that can be expressed in terms of the advanced Green’s functions of the system, 


Gcps = 27172 \G eii h 2 1 2 (18) 

Glar-i = ^ii\Ge u hi | 2 (19) 

Gqp-i = 27 Z Irn{^ m (0)}|G ez ,e m | ( 20 ) 

Gec = 0 ( 21 ) 


where the arguments of the Green’s function and hi refer to the electron and hole positions of the site i. The 
different contributions to the conductance are Cooper pair splitting, Gcps , local Andreev and the quasiparticle 
processes through QD Z, Glar-i and G qp -i, and elastic cotunneling, Gec- 


NON-INTERACTING CASE 

In this section, we derive an analytical expression for the total conductance through a Cooper pair splitter for the 
non-interacting Hamiltonian (with Ui = 0 in Eq. 0). Using standard techniques |2j we construct the exact advanced 
Green’s function from the Hamiltonians introduced in the previous section. For this purpose it is convenient to define 
the Nambu basis 


^ (t) d\^ d)^ d 2 ^ dii dm\, ^ 24 , 


jt A 


-1 


( 22 ) 


where £ t’ stands for the Heisenberg picture used for the dot operators. Then, we obtain the Green function 

-1 


g m 


(uj - e 1 - * 7i -tim 0 0 0 0 \ 

^1 m UJ 9mA) O m 0 ^indiA) 0 

0 — t 2m UJ — € 2 - *72 0 0 0 

0 0 0 uj + e i — * 7i / i ,„ 0 

0 Ai n d(uj) 0 ^1 m UJ T ^m 9m(uj) ^ 2 m 

0 0 0 0 t 2m w + e 2 - * 72 / 


(23) 


One can obtain analytical expressions of the conductances from Eq. (23). However, these expressions are cumbersome 
and it is more instructive to expand Eq. (23) perturbatively in terms of the induced gap A, n ,/, that is, e m , ti m A, ru |. 
Thus, to first order in A in( i we obtain 


G = 


uj T H e 


u> + H h 


{Cj + H e )~ l -A + HeC^A + Hh)- 1 

-A + HhA'AA + He)- 1 A + Hh)- 1 


( 24 ) 


where we defined the matrix cu = uj JT . =1 5ij, and 


/ zfei - i'ji 

“F hm 

° \ 

/ 

H e /h = I T^lm 

=F £m -9m(u) 

T^2m 1 

and A = 

V 0 

T^2 m 

T<7> - *72/ 

\ 


0 0 0 
0 A ind 0 
0 0 0 
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Then, we obtain the corresponding conductances 

(25) 

and (26) 

(27) 

The j-index stands for the opposite dot of QD l. From these expressions we find that CPS has a Lorentzian and 
the LAR and qp contributions a Fan-type shape. From the LAR and qp expressions one finds that the asymmetry 
orientation depends on the sign of e m , which one can thus relate to the shape-parameter in a Fano-type resonance [3]. 
We note that Glar-i and G qp -i exhibit the same gate dependence except for a factor 2 in the interference term. This 
factor accounts for the two particles in the local Andreev reflection. This factor is the origin of the more asymmetric 
interference patterns in the LAR contribution of the deeper minima compared to the quasiparticle contribution. 


Gcps ~ 27172 


Glar-i ~ 


^lra^2m'Tr? 


4(^+7?) 


2\2 


G 


qp-,1 


2 ^m7;Im{ff m (0)} 

4(^+7) 


2\ ’ 


e m( e l + 7l)( e 2 + 72 ) 

tlm^m ( 


-t 2 

L jm 


1 + 


e m (e]+7]) / 

et 2 \ 
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